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We search for the flavor-changing neutral current decays D° — > fi + fj,~ and D° — > e + e~ , and for the 
lepton-flavor violating decays D° — > e fj/^ using 660 fb _1 of data collected with the Belle detector 
at the KEKB asymmetric-energy e + e~ collider. We find no evidence for any of these decays. We 
obtain significantly improved upper limits on the branching fractions: B (D° — > /i + < 1.4 x 
10" 7 , B (D° -> e+e-) < 7.9 x 10" 8 and B (D° ->■ e+fi') + B (D° -> /x+e") < 2.6 x KT 7 at 90% 
confidence level. 

PACS numbers: 13.20.Fc,11.30.Hv,12.15.Mm,12.60.-i 



The flavor-changing neutral current (FCNC) decays 
D° — > e + e~ and D° — > n + /J.~ [I] are highly suppressed 
in the standard model (SM) by the Glashow-Iliopoulos- 
Maiani mechanism [2J. With the inclusion of long dis- 
tance contributions the branching fractions can reach val- 
ues of around 10 -13 [3]. The SM short distance Feynman 
diagrams for the D a — > decay are shown in Fig. [l] 

The lepton-flavor violating (LFV) decays D° — > e ± \i^ are 
forbidden in the SM, but are possible in extensions of the 
SM with nondegenerate neutrinos and nonzero neutrino 
mixings and are expected to be of the order of 10 -14 [3] 
in these scenarios. All these predictions are orders of 
magnitude below the current experimental sensitivity. 




FIG. 1. The SM short distance Feynman diagrams for the 
D° — > )i + fx~ decay. 

In certain new physics (NP) scenarios, FCNC branch- 
ing fractions can be enhanced by many orders of mag- 
nitude. For example, i?-parity violating supersymmetry 
can increase the branching fractions of D° — » e + e~ and 
D° u p to 10~ 12 and 1CT 8 , respectively [I]. 

The latter prediction is close to the current experimen- 
tal sensitivity. As another example, so far unobserved 
leptoquarks were suggested as a possible explanation of 
the small discrepancy between the measured value of the 
D s meson decay constant and the prediction of lattice 
QCD [5]. Leptoquarks could also enhance the D° —> 
£ + £~~ branching fraction. In order to explain the mea- 
sured Dj — » [i^v width by a leptoquark contribution, 
and comply with other constraints arising from charm 
meson decays, B (D — > should be enhanced to 

8 x 10 -7 [5J. The above examples demonstrate the im- 
portance of FCNC and LFV decays searches in the explo- 
ration of possible NP contributions. It should be noted 
that charm FCNC and LFV decays probe the couplings 



of the up-quark sector in contrast to B or K meson de- 
cays. 

In this paper, we report on a search for the decays 
D° fi + fi~, D° e+e~ and L>° -> e ± ^ using 
660 fb _1 of data recorded in e + e~ collisions at the center- 
of-mass (CM) energy of the T{AS) resonance and 60 MeV 
below by the Belle detector at the KEKB collider. 

The Belle detector, which is described in detail else- 
where [7J, is a large-solid-angle magnetic spectrometer 
that consists of a silicon vertex detector (SVD) [8], a 50- 
layer central drift chamber (CDC), an array of aerogel 
threshold Cherenkov counters (ACC), a barrel- like ar- 
rangement of time-of- flight scintillation counters (TOF) , 
and an electromagnetic calorimeter composed of CsI(Tl) 
crystals (ECL) located inside a superconducting solenoid 
coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside the coil is instrumented to detect 
if£ mesons and to identify muons (KLM). Two inner 
detector configurations were used. A beam pipe with 
a radius of 2.0 cm and a 3- layer silicon vertex detector 
were used for the first sample of 155 fb~ , while a 1.5 cm 
beampipe, a 4-layer silicon detector, and a small-cell in- 
ner drift chamber were used to record the remaining data 
sample. 

In this measurement only D° mesons coming from c - 
quark production in the continuum e + e~ — > cc process 
are considered. The inclusion of D° mesons from B de- 
cays would result in a higher combinatorial background. 
We normalize the sensitivity of our search to topologi- 
cally similar D° — > h + tt~ decays; this cancels various 
systematic uncertainties. The D° — > £ + £~ (£ = e or /i) 
branching fraction is determined by 

B (D° -> £+£-) = N te f (1) 

where Nju is the number of reconstructed D° — > £ + £~ 
decays and / is defined as 

/= J- 7T+7T-) (2) 

Here B (D° -> tt+tt") = (1.397±0.027) x 10~ 3 is the well- 
measured D° — > tt + it~ branching fraction 9 , is the 
number of reconstructed D° — > 7r + 7r~ decays, and eu and 
e n7r are the reconstruction efficiencies for D° — > £ + £~ and 
D° —> 7r + 7r~ decays, respectively. 
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First, a general event selection is performed that is, 
apart from the particle identification criteria, the same 
for all data samples. Later in the analysis, tighter opti- 
mized criteria specific for each decay mode are used. 

We use D° mesons from the decay D* + —> D°ir+ with 
a characteristic low momentum pion, since this consider- 
ably improves the purity of the D° — » £ + £~ and tt + it~ 
samples. Each charged track is required to have at least 
two associated vertex detector hits in each of the two 
measurement coordinates. To select pion and lepton can- 
didates, we impose standard particle identification crite- 
ria. Charged pions are identified using dE/dx measure- 
ment from the CDC, Cherenkov light yields in the ACC, 
and timing information from the TOF |10) . Muon identi- 
fication is based on the matching quality and penetration 
depth of associated hits in the KLM Electron iden- 
tification is determined using the ratio of the energy de- 
posit in the ECL to the momentum measured in the SVD 
and CDC, the shower shape in the ECL, the matching 
between the position of charged track trajectory and the 
cluster position in the ECL, the hit information from the 
ACC and the dE/dx information in the CDC [H]. The 
muon and electron identification efficiencies are around 
90% with less than 1.5% and 0.3% pion misidentification, 
respectively, whereas the pion identification efficiency is 
around 83%. D° daughter tracks are refitted to a com- 
mon vertex, and the D° production vertex is found by 
constraining the D° trajectory and the ir s track to origi- 
nate from the e + e~ interaction region; confidence levels 
exceeding 10~ 3 are required for both fits. A D* + mo- 
mentum greater than 2.5 GeV/c in the CM frame of the 
collisions is required to reject D-mesons produced in B- 
meson decays and to suppress combinatorial background. 

Candidate D° mesons are selected using two kinematic 
observables: the invariant mass of the D° decay prod- 
ucts, M, and the energy released in the D* + decay, 
q = (M D ,+ — M — m^c 2 , where M D ,+ is the invari- 
ant mass of the D°ir s combination and is the ir + 
mass (SJ. We require 1.81 GeV/c 2 < M < 1.91GeV/c 2 
and q < 20 MeV. 

According to Monte Carlo (MC) simulation based 
on EvtGen [13] and GEANT3 QJ], the background 
in D° — > £ + £~ decays originates predominantly from 
semileptonic B decays (80%) and from D° decays (10%). 
The background events can be grouped into two cate- 
gories based on their M distribution: (1) a smooth com- 
binatorial background, and (2) a peaking background 
from the misidentification of D° —> ir + ir~ decays. The 
decay D° — > ir + ir~ contributes to the background when 
both pions are misidentified as muons or as a muon- 
electron combination. MC studies showed that misidenti- 
fication of a K meson as a lepton does not produce a peak 
inside the mass region 1.81 GeV/c 2 < M < 1.91 GeV/c 2 . 

The signal efficiencies e« and e„ are evaluated using 
signal MC simulation, which is also based on EvtGen 
and GEANT3, but includes final-state radiative effects 
(FSR) simulated by PHOTOS [15;. Since we find dif- 
ferences between the widths of the D° —> 7r + 7r~ signal 



in the data and MC simulation, we perform fits to the 
M and q distributions to obtain scaling factors for the 
signal widths, and then tune the shapes in the MC sim- 
ulation by correcting M and q for each MC signal event 
by M' = to + (M - m )f m and q' = q + (q - q )f q . 
Here, too and qo denote the nominal D° mass and the 
nominal energy released in the D* + decay, respectively, 
and f m = 1.17 and f q = 1.28 are the corresponding scal- 
ing factors. Another difference, a small shift of —0.2 GeV 
in the mean of the distribution of the missing energy of 
the event (E m i ss ) between data and MC simulation is ob- 
served; we correct the MC distribution by subtracting, 
for every signal event, the above value from its -E m i ss - 
We construct E ni i SS from the difference between the beam 
energy and the sum of the energies of all four vectors of 
photons and charged tracks, which are assumed to be pi- 
ons. The constants derived from D° — > 7r + 7r~ are used to 
correct D —> £ + £r MC events. The uncertainties of the 
tuning procedure are included in the systematic error. 

In order to avoid biases, a blind analysis technique has 
been adopted. All events inside the D° signal region of 
| AM | < 20MeV/c 2 and \ Aq\ < 1 MeV were blinded until 
the final event selection criteria were established. Since 
D° — > £ + £~ decays are not expected to be observed at 
the current sensitivity, we optimize the selection criteria 
to obtain the best upper limits; we maximize the figure 
of merit, T — ^u/N\j\^, where eu is the efficiency for 
detecting D° —> £ + £~~ decays obtained from the tuned 
signal MC simulation and -/Vtjl is the Poisson average of 
Feldman-Cousins 90% confidence level upper limits on 
the number of observed signal events that would be ob- 
tained with the expected background and no signal |16j . 
The average upper limits Ajjl are calculated from the 
number of generic MC background events, surviving the 
selection criteria and scaled to the data size. The sample 
corresponds to 6-times the statistics of the data. 

For the optimization we select the following variables: 
signal region size (AM, Aq), -B m i SS , and minimal lep- 
ton identification probabilities. The quantities AM and 
Ag are measured relative to the nominal D° mass and 
nominal energy released in the D* + decay, respectively. 
The signal region in M is allowed to be asymmetric 
with respect to the nominal D° mass; for the \x\x de- 
cay mode this provides some suppression of misidenti- 
fied D° — > 7r + 7r~ decays, since their invariant mass dis- 
tribution peaks about 2 standard deviations below the 
D° mass; for the ee and e/x modes an asymmetric re- 
quirement accounts for the low mass tail due to electron 
bremsstrahlung. The requirement on the maximal al- 
lowed missing energy in the event is chosen to suppress 
background from semileptonic B decays; these events 
have larger missing energy due to undetected neutrinos. 
We found a broad maximum in T for the lepton iden- 
tification probability, hence we repeated the procedure 
at fixed lepton identification criteria, optimizing only the 
size of the signal region and the maximal allowed miss- 
ing energy in an event. The results are summarized in 
Table |H 
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TABLE I. Optimal selection criteria. The requirements on 
M are asymmetric and are given as lower and upper bounds 
on AM. 



I / / / Signal at 90% CL. gtffigi D°—*n + x~ Combinatorial 



Mode 


AM 


Aq 


-f^miss 




[MeV/c 2 ] 


[MeV] 


[GeV] 




(-8,19) 


±0.48 


1.4 


e + e~ 


(-27, 14) 


±0.40 


1.0 




(-13,15) 


±0.46 


1.0 



To estimate the number of combinatorial background 
events in the signal region, the sideband region \Aq\ > 
1 MeV is used. This region is chosen to reduce the sta- 
tistical error and to exclude possible signal events and 
misidentification from D° — > ir + ir~ decays. The com- 
parison of data and MC simulation shows good agree- 
ment in the combinatorial background distribution in this 
region. The distribution is parametrized as f(M,q) = 
A(l — BM)^/q, where the parameters A and B are deter- 
mined from a fit to the generic MC sample. The number 
of combinatorial background events in the signal region 
is calculated as N^f° = p x N sidc , where iV sidc is the 
number of events found in the sideband region and p is 
the expected ratio of events in the signal and sideband 
region determined by integration of f(M, q). 

The peaking background in the signal region due to 
misidentification of D° — s- it + tt~ decays is estimated from 
the reconstructed D a —> n + ir~ decays found in data by 
replacing the pion mass with the lepton mass and by 
weighting each event by 



u(pi, cos 9i)u(p 2 , cos 6> 2 ) 
V (pi , COS 6>i ) V (p 2 , cos 6» 2 ) 



(3) 



where px t % and #1,2 are the momenta and polar an- 
gles of the outgoing pions and where u and v are the 
pion-lepton misidentification probability and pion iden- 
tification efficiency, respectively. The misidentification 
probabilities and efficiencies are measured in data using 
D* + — » D° 7r+, D° — > K~ir + decays, binned in particle 
momentum p and cosine of polar angle. 

The estimates for the number of background events 
in the signal region are summarized in Table [TTJ. The 
misidentification of D° — > it + tt~ contributes significantly 
only to the D° — > n + ijT decay channel (1.8 events). The 
uncertainties in the background estimates listed in Ta- 
ble |n] include the statistical error on the number of side- 
band region events and the uncertainties in the combina- 
torial background parametrization, while the uncertainty 
in the peaking background estimation is negligible. 

The invariant mass distributions after applying the op- 
timized event selection criteria are shown in Fig. [2j In the 
signal region we find two candidates in the D° — > fi + /i~ , 
zero candidates in the D° — > e + e~, and three candidates 
in the D° — > e ± fi T decay mode; the yields are consistent 
with the estimated background. 

A binned maximum likelihood fit is used to deter- 
mine the yield of D° —> ir + ir~ candidates for the nor- 
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FIG. 2. The dilepton invariant mass distributions for (a) 
D° -> fi+fi-, (b) D° -> e+e", and (c) D° -> e ± ^ . The 
dashed vertical lines indicate the optimized signal window. 
Superimposed on the data (open histograms) are the esti- 
mated distribution for combinatorial background (filled his- 
togram) , the misidentification of D° — > iy + tt~ (cross-hatched 
histogram), and the signal if the branching fractions were 
equal to the 90% confidence level upper limit (single hatched 
histogram) . 



malization. We fit the invariant mass distribution using 
the same kinematic selections as for individual leptonic 
modes, except for the criteria on M. The fit function 
is the sum of two Gaussian distributions with the same 
mean and an FSR tail for the signal, and a first-order 
polynomial for the background. The shape of the FSR 
tail and its relative normalization are taken from the cor- 
responding signal MC simulation. The number of recon- 
structed D° mesons in the ir + ir~ mode is found to be 
51.2 x 10 3 , 44.1 x 10 3 , and 46.0 x 10 3 , using selection 
criteria for jijj,, ee, and efj, modes, respectively. The in- 
variant mass distribution of D° — > 7r + 7r~ using the /x/i se- 
lection criteria with the fit curve superimposed is shown 
in Fig. [3] The relative uncertainties on N^n are around 
0.5%. 

The signal efficiencies are determined from the tuned 
signal MC simulation. In addition, event weighting is ap- 
plied to compensate for small differences in lepton and 
pion identification efficiencies between data and MC sim- 
ulation. The correction factors for lepton identification 
were obtained using 77 — > £ + £~ and B — > XJ/ip(-^ 
£ + £~) decays. The signal efficiencies are found to be 
between 5% and 7% for £ + £~ decays and about 11% for 
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FIG. 3. The invariant mass distribution of D — > n n with 
the fit superimposed using [i/i selection criteria. 



TABLE II. Summary of the number of expected background 
events (Nbk a ), number of observed events (N) in the sig- 
nal region, the reconstruction efficiencies (e« and £-„-•„-) of the 
D° -> £ + r and D° — > tv + tt decays, the factors / and the 
branching fraction upper limits at the 90% confidence level. 





D" ->■ (j+fj,- 


D" -> e+e" 




Nbkg 


3.1 ±0.1 


1.7 ±0.2 


2.6 ±0.2 


N 


2 





3 


e«[%] 


7.02 ±0.34 


5.27 ±0.32 


6.24 ±0.27 


e™[%] 


12.42 ±0.10 


10.74 ±0.09 


11.22 ±0.09 


/[icr 8 ] 


4.84(1 ± 5.3%) 


6.47(1 ± 6.4%) 


5.48(1 ± 4.8%) 


UL [10-' 


] 1-4 


0.79 


2.6 



7r + 7r~ decays. The uncertainties in eu are estimated to 
be 0.3% and include contributions from MC statistics 
(0.2%), lepton identification efficiency corrections (0.2%), 
and MC tuning (0.1%). The uncertainty in e T7r is smaller 
(0.1%), because of a larger MC sample, better known 
pion efficiency corrections and a negligible contribution 
from MC tuning, since a wider range in M is used. 

From the number of reconstructed D° — > ir + ir~ de- 
cays, from the efficiency ratio, and from the known 
D° —> 7r + 7r~ branching fraction the factors / are calcu- 
lated with Eq. [2l The relative uncertainties are around 
5% (see Table In]) and include the errors on N n7r , e«, 
e-xTr and the D° — > 7r + 7r~ branching fraction, summed in 
quadrature. 

Finally, the branching fraction upper limits (UL) are 
calculated using the program pole . f [17], which extends 
the Feldman-Cousins method |16j by the inclusion of sys- 
tematic uncertainties. We find that the inclusion of sys- 
tematic uncertainties produces nearly the same result as 
the standard Feldman-Cousins method. The results are 
summarized in Table In] Note that B (D° e ± fi T ) de- 



notes the sum B (D° ->■ e+ fx~) + B (D° p+er). 

In summary, we have searched for the FCNC decays 
D° — > /i + pT and D° — > e+e~, and the LFV decays 
D° — > e ± /j, T using the Belle detector and have found 
no evidence of these decays. The upper limits on the 
branching fractions at the 90% confidence level are 

B(D° < 1-4 x 10~ 7 , 

B(D° ->■ e + e-) < 7.9 x 10~ 8 , 
B(D° -> e ± f i T ) < 2.6 x 10~ 7 . 

Previously, the best upper limits on these decays were 
published by the BaBar Collaboration [TB] using 122 fb _1 
of data. Our results improve these limits by a factor of 9 
for D° —> /i + /x _ decay, by a factor of 15 for D° —> e + e~ 
decay, and by a factor of 3 for D° —> e ± ^i =F decay. Re- 
cently, the CDF collaboration reported a preliminary 
result on the UL for the D Q — » branching frac- 

tion [19j : our result is lower by a factor of 3 and can 
further constrain the size of certain i?-parity violating 
couplings. It also strongly disfavors a leptoquark con- 
tribution [5] as the explanation for the anomaly in the 
measured Df —> p + v width [20 . 
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Throughout this paper charge-conjugate modes are in- 
cluded. 
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